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SIMULATION  OF  ELECTRON  BEAM  TRANSPORT  IN 
ION-FOCUSED  REGIME  CONDITIONING  CELLS 


I.  INTRODUCTION 

The  resistive  hose  instability^'^^  is  the  primary  factor  limiting  propagation  of  high 
current  relativistic  electron  beams  through  dense  gas  in  a  self-pinched  mode.  The  hose 
instability  is  a  macroscopic  m  s  l  kink  mode  which  disrupts  the  beam  by  amplifying 
transverse  beam  displacements.  The  quantity  r  =  t  -  zjc,  which  measures  the  time  into  the 
beam  pulse  at  a  fixed  axial  position  z,  is  an  q>proxunate  constant  of  the  motion  for  these 
beams,  and  dms  most  theoretical  treatments  lutve  used  z  and  r  (or,  equivalently,  $  s  cr)  as 
the  independent  variables  instead  of  z  and  t.  Theory  and  simulation  indicate  that  die  mode  is 
convective  in  the  beam  frame,  so  that  at  any  given  time  r  into  the  beam  pulse,  the  transverse 
displacement  Y(z,r)  grows  to  a  maximum  value  Y,^(r)  at  some  propagation  distance 
z  s  z^(r),  and  then  decays  as  the  beam  propagates  to  larger  z. 

The  most  unstable  modes  have  wavenumber  k  ■<  0.7k^  =  0.7(2irA^).  where  k^  and 
are  the  betatron  wavenumber  and  wavelength,  defined  by^‘^ 
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X^  H  2ir/k^  S  2»aj^(7l^/I^)*'  ,  (1) 

where  7  is  the  usual  Lorentz  factor,  a,,  is  a  characteristic  beam  radius,  and  I^  is  the  net  current 
(beam  current  I,,  less  plasma  return  current  Ip).  The  quantity  e  mc^/e  is  17  kA, 
independent  of  beam  parameters.  (The  Alfven  limiting  cunent  I^  =  jiyl,,.)  If  I),,  I„,  energy 
E  =  ymc^,  and  a^,  are  all  constant  in  r,  then  k^(r)  is  constant  throughout  the  beam,  and  the 
instability  can  grow  coherently  from  the  head  to  the  tail  of  the  beam.  Theory  and  simulation 
then  predict  many  orders  of  magnitude  of  hose  growth  for  typical  cases  of  practical  interest. 
In  experiments,  "raw"  beams  injected  from  a  beam  source  into  full  density  air  are  almost 
invariably  disrupted  by  hose  instability  after  propagating  a  few  betatron  wavelengths. 

For  many  purposes,  the  evolution  of  the  beam  radius  a|,(z,r)  can  be  adequately 
described  using  an  envelope  model  such  as  the  Lee-Cooper  equation.  This  model  asstimes 
that  the  beam  expands  or  contracts  self-similariy  with  the  radius  a|,  given  by 
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where  is  the  normalized  emittance  and  U  =  <k|r^  =  AJyrax?  represKtnts  Av*;  pinching 
effect  of  self-magnetic  fields.  A  pinched  equililmum  exists  when  =  T^SbU.  Equation  (2) 
assumes  the  beam  to  be  ultrarelativistic  (7  »  1)  and  paraxial.  The  paraxial  approximation 
assumes  that  the  transverse  vdocity  Vj_  of  beam  particles  is  much  less  than  tite  axial  velocity 
Vj  «  c.  The  Lee-Cooper  model  also  formally  defines  a^,  to  be  the  rms  radius  and  e  ^  to  be  the 
associated  rms  emittance.  In  addition,  if  the  net  current  and  beam  current  densities  have 
different  radial  profiles,  should  be  replaced  by  die  "effective”  current  a  radially- 

weighted  measure  of  the  net  current  within  the  beam. 

If  the  beam  characteristics  are  such  that  increases  monotonically  as  a  function  of  r, 
the  coherence  of  the  hose  mode  can  be  destroyed,  so  that  amplification  of  the  mode  is 
substantially  reduced.^  It  has  been  recognised  for  some  time  that  this  can  be  accomplished 
by  tapering  the  propagating  beam  fiom  a  broad  radius  at  the  head  to  a  much  narrower  radius 
at  the  tail.^'*^  It  is  readily  apparent  from  the  pinched  equilibrium  condition  that  it  is 
necessary  to  "tailor"  the  emittance  at  injection  to  insure  that  a|,(r)  will  be  tapered  during 
propagation.  Injecting  an  initially  tapered  af,(T)  into  a  dense  gas,  without  emittance  tailoring, 
will  not  lead  to  a  Uqiered  X^(t)  profile  that  is  maintained  during  propagation. 

Both  analytic  models^**  and  simulations^’^’**^^^  clearly  demonstrate  the  efiectiveness  of 
emittance  tailoring,  provided  the  a|,(  r)  profile  is  optimized  so  as  to  suj^ress  hose  throughout 
the  beam,  and  not  just  in  the  head  or  just  in  some  odier  part  of  the  beam. 

Even  with  a  beam  that  is  properly  emittance  tailored,  several  e-folds  of  hose 
amplification  can  be  expected  in  typical  cases.  Thus,  it  is  also  necessary  to  limit  the 
amplitude  of  initial  transverse  perturbations,  present  in  the  beam  at  the  point  of  injection  into 
dense  gas,  to  a  small  fraction  of  the  beam  radius.  [Actually,  these  perturbations  can  be  either 
displacements  X^fr),  Y^Ct)  or  transverse  velocities  (referred  to  as  "drifts")  dX^Crj/dz, 
dYg(r)/dz.]  The  frequency  content  of  these  perturbations  is  also  important.  At  any  given 
point  in  the  beam,  perturbations  with  frequency  conparable  to  r),  where 
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r^(T)  3  ira^(r)a^(T) /2c^ 

is  the  dipole  magnetic  decay  length,  induce  strong  hose  growth.^  Hete  o^lr)  is  the  on-axis 
ctmductivity,  which  rises  during  die  eaiiy  part  of  the  pulse  and  later  saturates  due  to  elecmm 
attachment  or  recombination.  High  frequency  perturbations  in  the  beam  head,  where  is 
short,  are  particularly  dangerous.  Low  frequency  displacements  do  not  couple  well  to  hose 
instability  and  are  relatively  benign. 

In  recent  years,  conditioning  techniques  have  been  developed  to  tailor  the 
emittance^*’^*^^’^®'^^  and  center  the  beam  to  mduce  initial  transverse  displacements.^^ 
These  techniques  have  been  implemented  in  experiments  with  the  and  ATA^** 

accelerators  at  Lawrence  Livermore  National  Laboratory,  SuperlBEX®’^^*^^*^^*^®  and 
PULSERAD^*  at  the  Naval  Research  Laboratory,  and  IBEX  and  RADLAC  at  the  Sandia 
National  Laboratories.^^  These  experiments  offer  dramatic  evidence  of  the  effectiveness  of 
beam  conditioning  in  stabilizing  hose. 

The  most  widely  used  technique  for  emittance  tailoring  involves  the  lue  of  a  "passive 
ion-focused  regime  (IFR)"  conditioning  cell.*®'^^  The  ion-focused  regime  refers  to 
propagation  of  a  highly  relativistic  beam  al<mg  a  plasma  channel  whose  charge  density  is  less 
than  that  of  the  beam.^®"^^’^*'^^  Since  7  »  1,  the  rq>ulsive  force  exerted  on  the  beam  by  its 
own  radial  electrostatic  field  E,  is  essentially  canceled  by  the  self-pinch  force  exened  by  the 
beam’s  azimuthal  magnetic  field  B^.  The  plasma  electrons  are  ejected  radially  by  E^,  leaving 
behind  a  positive  ion  column  which  pinches  the  beam  electrostatically.  The  background  gas 
density  must  be  sufficiently  low  that  the  plasma  electrons  are  collisionless  and  can  be  easily 
ejected.  If  the  plasma  charuiel  is  produced  by  the  beam  itself  through  impact  ionization  of  an 
initially  neutral  gas,  the  process  is  referred  to  as  passive  IFR.  If  the  gas  density  and  beam 
pulse  length  are  in  an  appropriate  regime,^’®*“'^®*^^’^*'“  the  ion  channel  density  increases 
significantly  during  the  beam  pulse,  and  therefore  the  beam  radius  tapers  down  throughout 
the  pulse.  This  radius  taper  may  be  converted  to  an  emittance  taper  by  passing  the  beam 
through  a  thick  scattering  foil  on  exit  from  the  IFR  cell.^*^*^^^  This  process  is  die  primary 
subject  of  the  present  paper,  we  shall  present  a  variety  of  axisymmetric  simulation  studies 
that  demonstrate  the  dependence  of  die  process  on  a  number  of  beam  and  control  parameters. 
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IFR  propagation  has  been  studied  in  other  contexts  by  several  investigators,  using  analytic 
models, axisynunetric  simulations,^^'^®*^^  and  3-D  simulations^®'^^  to  study 
non-axisymmetiic  effects  such  as  ion  hose  instability  and  magnetic  erosion. 

A  variety  of  techniques  have  been  used  to  center  the  beam  and  thus  reduce  the 
perturbations  which  "seed"  the  instability.  Passive  IFR  cells  exert  a  centering  force  due  to 
wall  image  forces,  but  the  effect  is  weak  if  the  wall  radius  is  considerably  larger  than  the 
beam  radius,  as  is  necessary  if  the  cell  is  intended  to  induce  radius  tapering.  However,  it  is 
also  possible  to  propagate  a  beam  in  an  "active"  IFR  cell,  wherein  some  external  ionization 
source  (e.g.,  a  laser^^"^^  or  a  magnetically-guided  low  energy  electron  beam^^)  establishes  die 
plasma  channel,  and  conditions  are  such  that  the  beam  itself  does  not  further  ionize  the  gas 
significantly.  This  type  of  cell  exerts  a  strong  centering  effect  on  the  beam,  and  damps 
transverse  perturbations  through  phase  mixing,  provided  the  density  profile  n^fr)  of  the  ion 
channel  is  not  constant  within  the  beam.  A  third  technique  is  to  propagate  the  beam  along 
either  a  current-carrying  narrow  wire^^*^*  or  an  electrostatically-charged  resistive  wire.^^*^’ 
The  highly  anharmonic  focusing  force  due  to  the  wire  centers  the  beam  and  strongly  danqis 
transverse  perturbations,  but  can  lead  to  undesirable  emittance  growth  if  parameters  are  not 
chosen  properly.  Conditioning  cells  based  on  all  of  these  techniques  have  been  employed,  in 
tandem  with  passive-IFR  tailoring  cells,  in  recent  experiments.  Figures  1  and  2  illustrate  the 
configurations  that  have  been  employed.  Figure  la  shows  a  simple  passive-IFR  tailoring 
cell,  with  a  thin  entrance  foil  FI  and  a  thicker  exit  foil  F2  for  scattering.  Figure  lb  shows  this 
cell  with  an  upstream  magnetic  focusing  lens,  a  configuration  that  was  employed  in  the  ATA 
experiments  to  provide  additional  control  over  the  radius  taper  a|,(r)  induced  in  the  cell.^ 
Figure  2  shows  three  configurations  for  a  passive-IFR  tailoring  cell  in  tandem  with  a 
centering  cell:  a  laser-IFR  cell  preceding  a  passive-IFR  cell,  often  used  in  the  ATA 
experiments  (Fig.  2a);  a  nanow  high-density  passive-IFR  centering  cell  preceding  a  broad 
lower-density  passive-IFR  cell,  used  in  earlier  SuperlBEX  experiments (Fig.  2b);  and  a 
passive-IFR  tailoring  cell  followed  by  a  wire  cell,  used  in  more  recent  SuperlBEX 
experiments  (Fig.  2c).  We  shall  present  simulation  studies  of  all  of  these  configurations. 

It  should  be  noted  that  if  magnetic  focusing  is  used  throughout  the  beamline,  as  was 
the  case  in  some  ATA  experiments,  the  focusing  system  can  be  tuned  to  minimize  the  spatial 


4 


pertuibations  X^Cr),  Y„(t)  at  the  point  of  injection  to  the  propagation  chamber.  However, 
doing  so  does  not  effectively  si^ress  hose  growth,  since  the  procedure  does  not  minimize 
the  transverse  velocities  dX^/dz,  dY^/dz  at  the  same  locati<m. 

The  organization  of  the  paper  is  as  follows.  In  Sec.  U,  we  present  studies  of  passive- 
IFR  tailoring  cells,  alone  or  in  combination  with  focusing  lenses  or  foils.  In  Sec.  m,  we 
present  studies  of  tandem  tailoring/centering  cells.  Section  IV  discusses  the  use  of  passive- 
IFR  tailoring  cells  in  multiple  pulse  experiments  such  as  the  Multi-Pulse  Propagation 
Experiment  (MPPE)  on  ATA.^  The  conclusions  are  summarized  in  Sec.  V.  The  Appendix 
describes  modifications  to  the  FRIEZR  simulation  code^^  that  was  used  for  all  studies 
presented  here. 
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n.  PASSIVE  IFR  TAILORING  CELLS 


A.  General  Description  of  a  Passive  IFR  Tailoring  Cell 

A  passive  IFR  consists  of  a  chamber  filled  with  low-density  gas 

(typically  a  few  mtorr)  in  which  the  beam  is  propagated  for  a  distance  equal  to  several 
betatron  wavelengths.  The  channel  strength,  and  therefore  the  pinch  force,  increase 
monotonically  with  r,  the  time  into  the  beam  pulse.  The  very  front  of  the  beam  head  is 
uiq>inched,  and  thus  expands  at  a  constant  ram  characteristic  of  free  expansion,  da^dz  -  Ad, 
where  Ad  is  the  thermal  microdivergence  in  beam  electron  velocity  angles.  The  weakly 
pinched  beam  head  expands  rapidly,  but  the  beam  may  need  to  propagate  a  long  way  before 
the  beam  head  reaches  equilibrium  at  a  large  radius.  In  fact,  at  any  given  propagation 
distance  z,  there  will  always  be  a  weakly  pinched  portion  of  the  beam  head  which  is  still 
expanding,  and  for  the  typical  short  lengths  of  IFR  conditioning  cells,  this  expanding  region 
may  constitute  a  substantial  fraction  of  the  pulse.  The  propagation  distance  required  for 
equilibration  decreases  as  a  function  of  r.  The  tail  quickly  reaches  an  equilibrium  radius  that 
(if  entrance  conditions  are  chosen  optimally)  is  C(»iq>arable  to  the  initial  radius.  If  the  cell  is 
long  enough,  the  beam  evolves  toward  an  equilibrium  in  which  the  radius  decreases 
monotonically  as  a  function  of  r.  The  normalized  onittance,  =  ytH^Ad,  typically  will  not 
develop  this  monotonic  tiq^*r  within  the  IFR  cell;  in  many  cases  it  is  qjproximately  true  that 
A^  decreases  as  a|,  increases,  so  that  remains  roughly  constant.  However,  at  the  end  of  the 
cell  the  beam  is  passed  through  an  exit  foil  which  scatters  the  electrons  sufficiently  to 
increase  A0  to  a  value  which  is  roughly  uniform  over  the  duration  of  the  pulse.  The  radius 
taper  is  thus  effectively  converted  to  an  emittance  taper. 

A  simple,  qualitative  model  of  passive  IFR  conditioning  may  thus  be  based  on  the 
following  assumptions:  (1)  That  is  constant  within  the  cell.  (2)  That  ion  channel  motion 
may  be  neglected.  (3)  That  the  beam  expands  self-similarly,  and  that  the  beam  and  channel 
radial  profiles  are  similar.  (4)  That  a|,(z,  r)  is  equal  to  the  smaller  of  the  equilibrium  radius 
a^q(r)  and  the  free  expansion  radius  a^^^z).  (S)  That  the  beam  is  neither  converging  nor 
diverging  at  injection,  i.e.,  that  da^dz  =  0  at  z  s  0.  The  resulting  analytic  model  is  presented 
more  completely  and  in  a  somewhat  different  form  in  Ref.  20,  and  is  an  extension  of  a 
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steady-state  model  by  Briggs.^  It  provides  considerable  insight  and  is  useful  for  scaling, 
although  all  of  these  assun^ions  are  of  limited  validity.  In  subsequent  sections,  we  shall 
elaborate  on  the  limitations  of  these  assumptions  and  use  particle  simulations  to  calculate 
beam  evolutitm  mote  exactly. 

With  these  assunq)tions,  we  proceed  to  calculate  the  beam  radius  and  emittance.  We 
assume  the  beam  is  injected  into  the  IFR  cell  with  radius  a^,  thermal  microdiveigence 
and  normalized  emittance  s  The  density  of  ions  created  by  beam-impact 

ionizaticm  is  given  by^*^ 

hnJhT  -  (4) 

where  n),(r,r,z)  and  nj(r,r,z)  are  the  beam  and  ion  densities,  respectively,  and  Vj  is  the  gas 
ionization  rate.  In  air  at  pressure  P,  Vj  •<  P(torr)  nsec*^  In  the  body  of  the  beam,  dte  radius 
a^(  t)  is  given  by  the  Bennett  pinch  condition^®'^*^ 

1/2 

•q  0  o  b 

The  ratio  f(r)  s  n^fr  =  0,r)/ni,(r  s  0,  r)  is  used  to  characterize  die  strength  of  the  channel. 
Equation  (S)  is  the  equilibrium  solution  to  the  Lee-Cooper  envelope  equation  (2)  with  I^ 
replaced  by  fl),.  The  radius  z^{z)  s  a|,(r  =  0,  z)  at  the  beam  head  is  specified  by  free 
expansion. 

The  free  expansion  limit  can  be  obtained  from  Eq.  (2)  with  U  =0,  and  reduces  to  the  second 
expression  in  Eq.  (6)  when  a^,  »  a^.^^  If  is  defined  as  some  suitable  reference  point  in 
the  body  or  tail  of  the  beam,  the  ratio  a^/a^^frp)  may  be  referred  to  as  the  radius  tailoring 
ratio. 

In  die  region  of  the  beam  that  has  reached  equilibrium,  the  microdivergence  is  given  by 

Ad  -  (fl. /tI  (7) 

•q  bo 
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and  increases  monotonically  with  t.  The  beam  radius  taper  is  converted  to  the  desired 
emittance  taper  by  passing  the  beam  through  a  thick  scattering  foil  located  at  z  =  Zf.  The  foil 
adds  a  scattering  angle  to  the  beam,*’  resulting  in  a  thermal  microdivergoice  given  by 


{z  )  -  . 

I  8C  eq 


(8) 


The  emittance  C||(Zf)  at  the  beam  head  after  passing  duough  the  foil  is  thus 


«  A8  z, 
o  sc  f 


(9) 


Sunilaily,  die  emittance  ejiZf)  in  the  beam  tail  is 


1/2 

e 

Ae"  +  — - 

0 

“b 

sc  7I 

oj 

1/2 


(10) 


The  ratio  e„(Zf,  0)/e„(Zf,  Tp)  *  ftom  (9)  and  (10)  defines  an  emittance  tailoring  ratio. 
To  complete  the  calculation,  it  is  necessary  to  specify  the  neutralization  fraction  f(  t).  For  a 
beam  with  constant  current,^  f  =  v^rP.,  For  arbitrary  7(t)  and  €^{t)  profiles  at  z  =  0, 
f  can  be  estimated  from® 


fir) 


V  .  * 

i  o 

2  47  I^(r) 

ID 


fr 


t 

o 


(11) 


As  an  example,  consider  beam  parameters  typical  of  the  ATA  Multi-Pulse  Propagation 
Experiment^  (MPPE);  7  =  20,  ag  =  0.5  cm,  «q  =  0.5  rad-cm,  propagating  60  cm  in  10  mtorr 
of  air  followed  by  a  scattering  foil  with  =  0.09.  The  ionization  rate  at  this  density  is 
Vj  =  10^  sec"*.  From  (4)  and  (9),  the  radius  a|,(Zf,  0)  and  emittance  0)  in  the  beam  head 
at  Zf  =  60  cm  are  3.0  cm  and  5.4  rad-cm,  respectively.  If  the  beam  is  30  nsec  long,  the  tail 
radius  in  equilibrium  from  (8)  is  0.49  cm.  The  final  emittance  from  (9)  is  ej(Zf)  =  1.0  rad- 
cm,  resulting  in  a  ratio  =  5.4.  This  tailoring  ratio  is  in  the  range  considered  suitable 
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for  reducing  the  growdi  of  die  resistive  hose  inability  to  an  acceptable  level  when  the  beam 
is  subseqaendy  propagated  in  dense  gas. 

B.  Overview  of  FRIEZR  Simulation  Code 

The  simple  model  outlined  in  the  previous  section  neglects  a  number  of  effects  which 
may  be  treated  more  precisely  using  time>dependent  envelope  codes  or  full  particle 
simulation  codes  such  as  FRIESl..^^  First,  there  is  a  transition  region  widiin  the  beam  pulse 
where  the  beam  is  weakly  pinched  and  has  not  yet  reached  an  equilibrium.  In  many 
experimental  situations,  the  tailoring  ceU  is  so  short  that  a  large  fraction  of  the  beam  has  not 
yet  reached  equilibrium  at  the  end  of  the  cell,  so  die  analytical  model  cannot  be  relied  on  for 
making  detailed  predictions  of  emittance  tailoring  profiles.  Full  simulation  codes  such  as 
FRIE21R  include  a  number  of  physical  efifects  which  are  usually  omitted  from  simpler 
models.  For  example,  the  FRIEZR  code  includes  beam  energy  changes  arising  from  the 
inductive  electric  field  and  does  not  make  the  ultrarelativistic  and  paraxial  approximations  in 
treating  beam  dynamics.  Particle  simulations  also  allow  the  beam  cturent  density  profile  to 
evolve  self>consistently  and  thus  do  not  make  the  self-similar  assumptions  used  in  envelc  ps 
models.  In  addition,  FRIEZR  can  treat  more  conq>licated  beamline  configurations  and  multi¬ 
stage  conditioning  cells.  These  beamline  configurations  may  include  vacuum  transport 
regions  with  solenoidal  fields  or  discrete  focusing  lenses  as  well  as  centering  cells  employing 
preionized  plasma  channels. 

In  the  remainder  of  this  pq>er,  we  concentrate  on  simulation  studies  that  are  chosen  to 
model  two  recent  series  of  experiments  conducted  with  the  ATA  accelerator^  at  Lawrence 
Liveimore  National  Laboratory,  and  the  SuperlBEX  beam  generatoi*’’’^^’^®'^’^  at  the  Naval 
Research  Laboratory.  In  diis  particular  series  of  ATA  experiments,  known  as  the  Multi-Pulse 
Propagation  Experiment  (MPPE),  a  train  of  up  to  five  10  MeV  pulses,  separated  by  several 
msec,^  was  conditioned  and  propagated.  (Earlier  single  pulse  experiments  at  ATA  used 
energies  in  excess  of  40  MeV.)  SuperlBEX  is  a  single  pulse  device  which  produces  a  4.5 
MeV  beam  at  currents  which  may  exceed  20  kA.  Table  I  lists  the  bean,  and  conditioning  cell 
parameters  used  (excq>t  as  noted  specifically  in  the  text)  in  all  of  our  simulations  of  these 
experiments.  In  all  of  the  simulations,  the  gas  in  the  IFR  tailoring  cell  was  assumed  to  be  air. 
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(Several  oAer  gases  as  well  as  air  were  used  in  the  e^qjeriments.)  The  plasma  electrons  were 
assumed  in  all  cases  to  be  ejected  instantaneously  from  die  channel. 


TABLE  1.  Beam  and  IFR  cell  parameters  for  analytical  and  simulation  exanqiles. 


Parameter 

Symbol 

Units 

ATA 

SuperlBEX 

Beam  energy 

y 

- 

21 

10 

Beam  current 

ho 

kA 

6 

15 

Initial  half-current  beam  radius 

cr.. 

0.50 

0.60 

Initial  rms  beam  radius 

WO) 

cm 

0.60 

0.72 

Current  rise  length 

’’r 

nsec 

10 

10 

Pulse  length 

nsec 

40 

26.7 

Initial  rms  emittance 

s 

a 

II 

0 

rad-cm 

0.48 

1.08 

Initial  half-current  emittance 

^ho 

rad-cm 

0.40 

0.90 

IFR  cell  radius 

cm 

7.5 

2.5/5 

IFR  cell  pressure 

p 

mtorr 

2-50 

20/5-100 

The  primary  quantity  of  interest  from  the  simulations  is  the  beam  radius  a,,(T)  at  the 
IFR  cell  exit  foil.  Up  to  this  point,  we  have  been  rather  cavalier  in  discussing  the  "beam 
radius"  Since  pinched  beams  generally  have  rounded  beam  current  density  profiles  Ji,(r), 
there  are  potentially  many  ways  to  define  a  characteristic  beam  radius  and  emittance  in 
particle  simulation  codes  or  experimental  data.  The  rms  radius  a^^,  is  the  quantity  that 
appears  in  the  Lee-Cooper  definition  of  emittance.  However,  the  rms  radius  can  be 
dominated  by  a  small  number  of  beam  electrons  in  a  "halo"  at  large  radius.  We  generally 
find  the  half-current  radius  to  be  a  more  useful  indication  of  the  beam  width.  For  any 
specific  radial  profile,  ^  well-defined  value,  e.g.,  ~1.2  for  a  gaussian.  If  this 

ratio  becomes  large,  it  is  an  indication  of  change  in  the  radial  profile,  typically  with  the 
growth  of  radial  wings. 
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The  beam  in  FRIEZR  is  defined  by  a  half-cuiient  radius  at  injection  aj^(z  s  0)  =  a^, 
nominal  beam  eneigy  curtent  current  rise  time  pulse  lengdi  fp  and  initial 

mw  emittance  e^.  Beam  simulation  particles  are  loaded  with  a  uniform  distributitm  in  die  die 
axial  variable  t,  but  the  charge  on  the  simulation  particles  is  weighted  by  W(t)  = 

In  most  cases,  the  currem  is  assumed  to  have  die  form 

I.  (t)  »  I  tanh{T/T  ),  (12) 

DO  T 

Particles  are  distributed  in  the  radial  direction  with  an  assumed  Gaussian  current  density 
profile  of  the  form: 

Jj^(r)  «  — I  ejcp(r^/ag) ,  (13) 

where  a^  and  Ij,  may  vary  with  r.  For  a  Gaussian,  a^  “  l-2ajy2.  [FRIEZR  has  options  for  a 
number  of  Ib(T)  and  Jj,(r)  profiles  other  than  those  in  Eqs.  (12)  and  (13).]  The  initial 
normalized  transverse  velocity  Vj/c  for  individual  simulation  particles  is  chosen  from  a 
Maxwellian  distribution  with  a  thermal  microdivergence  ^0^  which  is  related  to  the  initial 
emittance  by 


e  /y  a  (0)  , 
o  o  rms 


(14) 


The  studies  presented  here  do  not  include  7(r)  variations  which  may  arise  from  experimental 
voltage  variations  or  space  charge  depression.  Space  charge  depression  effects  ate  relatively 
unimponant  for  IFR  transport  in  the  ATA  regime  but  can  be  significant  for  SuperlBEX.^ 

The  FRIEZR  code  uses  an  unconventicmal  coordinate  syston  which  moves  at  die  speed 
of  light  with  r,  t  and  t  =  t  -  z/c  as  the  independent  variables.  (Note  that  much  of  the  previous 
literature  on  FRIEZR  and  related  simulation  codes  has  used  5  h  cr  rather  than  r  as  the  axial 
coordinate  for  the  distance  from  the  beam  head.)  The  z  coordinate  does  not  appear  explicidy 
in  the  model  equatitms,  so  z-dependent  beamline  changes  must  be  transformed  to  the  (r,  t,  t) 
variables  using  z  « t  -  cr.  Data  such  as  the  beam  radius  or  emittance  are  usually  plotted  as 
functions  of  r  at  fixed  z-locadons  in  order  to  facilitate  conqiatisons  with  time>resolved  beam 
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diagnostics.  The  implications  of  this  coordinate  system  and  other  numerical  details  are 
discussed  in  Ref.  47.  The  ^ypendix  includes  additional  discussion  of  numerics  as  well  as  a 
descriptim  of  features  added  specifically  to  treat  conditimimg  cells.  These  features  include 
beam  impact  ionization,  channel  ion  motion,  scattering  and  focusing  thin  foils,  discrete 
focusing  lenses,  and  z-dependent  wall  radius,  preionized  channel  doisity  and/or  gas  pressure. 

FRIEZR  has  been  used  to  study  inductive  erosion  in  preionized  channels,^^  wakefield 
effects  in  laser  propagation  through  a  plasma,^^  beam  transport  in  the  ATA^*  and  RLA^ 
accelerators,  beam  propagation  in  passive  IFR  cells,^^  and  both  ballistic  and  pinched 
propagation  in  heavy  ion  fusion  target  chambers.^ 

C.  Non-Diverging  Beam  Injected  into  a  Pa^ve-IFR  Cell 

In  this  section,  we  consider  a  simple  passive-IFR  cell,  as  shown  in  Fig.  la.  We 
continue  to  assume,  for  the  present,  that  the  beam  is  neither  converging  nor  diverging  at 
injection,  i.e.,  that  da^dz  =  0  at  z  =  0.  In  so  doing,  we  are  neglecting  the  focusing  effect  of 
the  entrance  foil,  and  excluding  the  possibility  of  a  focusing  lens  iqystream.  We  consider  five 
cases  summarized  in  Table  n,  each  representing  a  specified  variation  of  the  ATA  beam  of 
Table  1.  The  beam  characteristics  at  the  output  end  of  the  tailoring  cell  (z  s  90  cm)  are 
summarized  in  Table  n  in  terms  of  a|^(r  s  0),  die  half-current  radius  of  the  beam  head;  and 
ai^(r  =  40nsec),  a^^gCr  =  40),  =  40),  the  half-radius,  rms  radius,  and  rms  normalized 

emittance  at  the  beam  tail. 

In  Case  1-A,  we  specify  that  the  channel  ions  ate  stationary,  and  the  initial  emittance  is 
roughly  matched  to  the  pinch  force  exerted  by  the  beam-ionized  channel.  The  tail  radius  thus 
remains  essentially  constant,  while  the  rest  of  the  beam  expands.  Figure  3  shows  that  the 
emittance  remains  reasonably  constant  throughout  the  beam.  Thus  the  assun^tions  of  the 
analytic  theory  of  Sec.  II.A  are  reasonably  well  satisfied.  Indeed,  Fig.  4  shows  that  the 
output  beam  radius  profile  -  90,  r)  agrees  well  with  the  analytic  theory,  particularly  in 
the  well-pinched  tail. 
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TABLE  n.  Sununaiy  of  FRIEZR  simulations  of  simple  passive  IFR  tailoring  cells.  Output 
quantities  (last  four  columns)  are  at  z  s  90  cm  and  r  -  0  (beam  head)  or  r  =  40  nsec  (beam 
tail). 


Case  Comment  P  a,y2(’'=0)  W^O)  e„(40) 


rad-cm 

mtorr 

cm 

cm 

cm 

rad-cm 

1-A 

Fixed  ions 

0.48 

5 

2.8 

0.48 

0.62 

0.50 

1-B 

LowP 

0.48 

5 

2.8 

0.38 

0.57 

0.55 

1-C 

HighP 

0.48 

20 

2.8 

0.31 

0.55 

1.05 

1-D 

High  «„,P 

0.96 

20 

4.2 

0.41 

0.61 

1.2 

1-E 

ao=  1.0  cm 

0.96 

20 

2.0 

0.57 

0.93 

1.95 

Case  1-B  is  similar  except  that  ion  dynamics  ate  included  (and  also  in  aU  the  remaining 
cases  discussed  in  this  piq>er).  In  the  beam  tail,  the  ions  collapse  toward  the  axis  under  die 
influence  of  the  beam’s  radial  electrostatic  field.^’^  For  f  «  1,  the  characteristic  time  scale 
for  ion  channel  coUtqise  is 

1/2  1/2 

T  *  (xa  /2c)  (m. /m  )  (I  /I  )  (15) 

c  o  X  e  o  D 

The  ion  collapse  increases  the  strength  of  the  pinch  force  on  the  beam,  and  thus  reduces  the 
beam  tail  radius  8^/2  at  z  s  90  cm  from  0.48  cm  in  Case  1-A  to  0.38  cm  in  Case  1-B.  This 
effect  is  generally  desirable,  although  ion  channel  motion  may  also  lead  to  beam  degradation 
through  ion  hose  instability Figure  5  shows  the  evolution  of  Siip_{T)  as  the  beam  is 
transported  through  the  IFR  cell  to  z  =  30, 90,  and  ISO  cm.  In  agreement  with  the  anal3mc 
model,  the  head  expands  freely  while  the  tail  rapidly  reaches  equilibrium.  There  is  a 
transition  region  (3  nsec  <  r  <  17  nsec)  where  the  beam  expands  at  a  reduced  rate.  Figure  6 
shows  the  expansion  rate  da^^dz  as  a  function  of  r  at  the  same  three  z-locations.  (In  this 
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particular  figure,  we  plot  the  rms  radius  rather  than  the  half-current  radius.)  At  r  =  0,  the 
beam  continues  to  e^and  at  a  constant  rate  equal  to  the  free  expansion  rate.  An  qrproximate 
equilibrium  is  reached  (dajg,,/dz  •<  0)  in  the  tail  of  the  beam,  with  the  equilibrium  region 
moving  forward  as  the  beam  propagates  further  in  z.  Thus  the  tapered  rms  radius  a^,(  r) 
continues  to  evolve  as  the  beam  propagates  through  a  cell  of  this  length.  The  "equilibrium" 
radius  in  the  beam  tail  actually  is  seen  to  expand  slowly.  This  occurs  for  two  reascms.  First, 
the  strength  of  the  channel  ionized  by  the  beam  gradually  drops  as  the  beam  head  expands, 
from  its  initial  value  fBVjr  for  a  pencil  beam,  toward  the  Briggs^  equilibrium  value  f«Vjr/2. 
Secondly,  the  collq)se  of  the  ion  channel  causes  the  beam  radial  density  profile  n|,(r)  to  foim 
a  mild  halo,  which  increases  the  ratio  from  its  initial  value,  and  thus  increases 

da^^^dz.  It  can  be  seen  in  Fig.  7  that  a^,/a|/2  in  the  tail  increases  to  1.46,  from  its  initial 
value  of  1.2. 

In  Case  1-C,  the  gas  density,  and  therefore  the  quantity  f,  which  is  an  index  of  the 
pinch  strength,  are  increased  by  a  factor  of  four.  The  simple  model  of  Sec.  IIA  would 
indicate  that  a^/j  is  reduced  by  a  factor  of  two.  However,  a  more  detailed  analysis^  reveals 
that  the  beam  cannot  collapse  to  an  equilibrium  radius  less  than  if  the  initial 

emittance  is  too  small  to  balance  the  pinch  force,  then  the  emittance  must  increase  as  the 
beam  equilibrates.  (The  opposite  is  not  true:  if  is  too  large  to  balance  the  pinch  force,  the 
beam  can  expand  with  little  or  no  emittance  increase.)  Indeed,  Fig.  8  and  Table  n  show  that 
Si^l2  collapses  to  0.62aQ,  and  Fig.  3  clearly  shows  the  emittance  increase  in  the  tail  that  is 
necessary  to  prevent  further  radial  collapse.  ITie  increase  in  a^g/aj^  in  the  tail,  from  its 
initial  value  of  1.2  for  a  gaussian  profile  to  its  final  value  of  1.8,  is  an  indication  of  the  profile 
deformation  (enhancement  of  the  radial  wings)  that  accompanies  the  emittance  increase.  The 
tailoring  induced  in  this  case  would  not  be  ideal  for  typical  propagation  applications,  since 
the  taper  is  restricted  to  a  small  region  r  ^  10  nsec  in  the  head,  leaving  the  remainder  of  the 
beam  essentially  untailoied. 

In  case  1-D,  the  initial  beam  parameters  are  similar  to  Case  1-C  except  that  the 
emittance  is  doubled.  As  a  result,  the  tail  is  initially  close  to  equilibrium,  and  the 
evolution  is  similar  to  Case  1-B  but  at  twice  the  radius.  However,  in  this  case  the  outer  pan 
of  the  beam  head  contacts  the  cell  wall  at  r  =  7.5  cm,  and  is  lost.  As  a  result,  decreases  in 
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the  head,  as  is  seen  in  Fig.  3.  Loss  of  pait  of  die  beam  head  limits  emittance  tailoring  and  is 
to  be  avoided  in  practical  inqilementations. 

Case  l-E  has  die  same  initial  emittance  as  Case  1-D,  but  the  initial  radius  a^  is  doubled 
while  the  initial  velocity  spread  is  halved.  As  in  Case  1-C,  the  initial  emittance  is 
insufficient  to  support  the  pinch  force  in  the  tail,  so  that  ajy2  collqises  in  the  tail  to  0.57ao, 
while  increases  to  1.7  (indicative  of  halo  formatim),  and  increases  in  the  tail,  as 

is  evident  in  Fig.  3. 

D.  Effects  of  IFR  Entrance  Foil  and  Upstream  Focusing  Lenses 

In  previous  sections,  it  was  assumed  that  the  beam  was  neither  converging  nor 
diverging  at  entiy  to  the  IFR  cell,  i.e.  da^dz  s  0  at  z  =  0  for  all  beam  slices.  If  a  conducting 
foil  is  present  at  entiy  to  the  cell,  it  acts  as  a  focusing  lens  by  shorting  out  die  beam’s  space 
charge  field.^^’^^*^*  The  focusing  angle  decreases  with  7  and  is  roughly  proportional  to 
^(r,  r),  the  current  enclosed  within  radius  r.  The  focusing  strength  thus  varies  with  time  imo 
the  pulse  and  is  nonlinear  in  r.  Foil  focusing  thus  is  not  easily  controlled  and  tends  to  be 
disruptive.  However,  it  is  also  possible  to  introduce  one  or  more  magnetic  lenses  in  the 
vacuum  transport  region  ahead  of  die  IFR  entrance  foil.  These  lenses  can  be  precisely 
controlled,  and  can  be  used  to  partially  compensate  for  foil  effects  or  otherwise  to  fine>tune 
the  radius  profile.  This  method  was  extensively  and  successfully  employed  in  the  MPPE 
experiments.^ 

In  this  section,  we  present  simulations  of  the  conditioning  configuration  shown  in  Fig. 
lb:  a  focusing  lens  L  located  at  z  =  6  cm,  followed  by  a  vacuum  transport  regiirn,  and  then 
by  a  passive-IFR  tailoring  cell  with  an  entry  foil  FI  at  z  =  126  cm  and  a  thicker  exit  foil  F2  at 
z  =  176  cm.  The  beam  parameters  are  representative  of  the  ATA  experiments;  they  are  as 
listed  in  Table  I  except  that  in  this  section  the  beam  rise  time  is  set  to  r,  =  20  nsec,  and  (in  all 
cases  except  2-B)  =  0.24  rad-cm.  The  beam  is  injected  with  a  relatively  large  radius  aj„j. 

The  values  of  a^^j  and  of  the  lens  focal  length  f^^  are  specified  in  Table  in.  In  all  cases  excqn 
2-G  and  2-J,  the  entry  foil  FI  is  effectively  a  10  mil  carbon  foil,  scatteriitg  the  beam  through 
an  angular  spread  ~  0.045  radians. 
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The  results  of  the  simulations  are  shown  in  Table  m  and  Figs.  9-13.  The  effectiveness 
of  the  tailoring  process  may  be  estimated  by  comparing  the  radii  ay2  at  the  end  of  the  IFR 
cell,  given  in  Table  m  for  r  =  0,  20  and  40  nsec,  the  head,  middle  and  tail  of  the  beam, 
respectively. 

TABLE  m.  Summary  of  FRIEZR  simulations  using  an  upstream  focusing  lens  and  a  SO  cm 
long  passive  IFR  cell  widi  air  at  6  mtorr.  All  lengths  are  in  cm. 


Case 

k 

cm 

Vi 

cm 

cm 

av2(20) 

cm 

cm 

Comments 

2-A 

114 

1.5 

2.25 

1.48 

1.10 

Small  ajnj;  large  4 

2-B 

114 

1.5 

2.6 

1.54 

1.24 

Cq  =  0.48  rad-cm 

2-C 

114 

2.0 

2.25 

1.47 

0.84 

Small  ajnj;  large 

2-D 

114 

2.5 

2.4 

1.59 

0.92 

f^  too  large 

2-E 

114 

3.5 

2.55 

2.05 

1.61 

Large  aj^j  and  fj. 

2-F 

100 

2.5 

2.65 

1.51 

0.64 

Near  optimum  f|^ 

2-G 

90 

2.5 

2.5 

0.99 

0.39 

=  0.031 

2-H 

90 

2.5 

2.95 

1.44 

0.57 

Near  optimum  f^ 

2-J 

90 

2.5 

3.4 

2.25 

1.05 

=  0.063 

2-K 

80 

2.5 

3.2 

1.51 

0.67 

Near  optimum  f^ 

2-L 

70 

2.5 

3.65 

1.80 

0.92 

f too  ^ort 

For  these  conditions,  focal  lengths  f^  >  100  cm  (Cases  2-A  through  2-E)  lead  to  poor 
tailoring.  Although  the  distance  between  the  lens  and  the  entrance  foil  is  120  cm.  Figure  9, 
which  plots  a^^(z)  for  three  beam  slices  taken  from  simulation  Case  2-D,  shows  that  beam 
space  charge  repulsion  extends  the  lens  focal  length  considerably,  particularly  in  the  high- 
current  part  of  the  beam  (r  =  20  and  40  nsec).  Because  the  entire  beam  has  a  large  radius  and 
is  converging  at  the  IFR  entry  foil,  the  beam  subsequently  overfocuses  within  the  IFR  cell 
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and  rebounds,  leading  to  a  rather  large  tail  radius  (^0.92  cm)  at  the  exit  foil,  and  insufficient 
titering.  Note  the  effect  of  foil  focusing  by  the  IFR  entrance  foil,  which  is  evident  in  Fig.  9 
as  a  sudden  change  in  the  slope  of  the  a^^(z)  trajectory  at  z  =  126  cm.  This  foil  heats  the 
beam  due  to  both  scattering  and  anharmonic  focusing.  Beyond  this  point,  the  beam  is 
pinched  by  the  ion  channel  whidi  slows  the  expansion  in  the  beam  body. 

Figure  10  plots  timiz)  for  the  same  three  beam  slices  when  the  focal  length  is  reduced 
to  90  cm  (Case  2-H).  The  t -A  nsec  and  20  nsec  beam  slices  now  focus  well  ahead  of  the 
entrance  foil  while  the  rs  40  nsec  beam  slice  focuses  just  ahead  of  the  foil.  This  variation  in 
the  focal  point  serves  to  increase  the  beam  head  radius  and  decrease  the  minimum  radius, 
thus  substantially  improving  the  radius  tzpct  after  transiting  the  tailoring  cell. 

Figure  1 1  plots  tLi/2(‘r)  just  after  the  final  scattering  foil  for  f^  =  70,  90  and  1 14  cm 
(Cases  2-L,  2-H  and  2-D).  The  90  cm  focal  length  produces  the  smaUest  beam  radius  in  the 
body,  while  the  other  two  cases  produce  nearly  identical  results  for  r  >  23  nsec.  The 
corresponding  final  eminance  profiles  are  shown  in  Fig.  12;  the  f^  =  90  cm  case  produces  an 
emittance  variation  of  a  factor  of  3. 

The  tapered  shq>e  of  the  beam  radius  a|/2(v)  is  often  quite  sensitive  to  the  thickness  of 
the  IFR  entrance  foil  FI.  This  effect  has  been  noted  in  IFR  conditioning  cell  experiments 
with  both  MPPE^  and  SuperlBEX.^^*^  Figure  13  plots  the  final  radius  tn/ziT)  for  Cases  2-G, 
2-H  and  2-K,  which  are  similar  except  for  thickness  of  the  IFR  cell  entrance  foil.  In  these 
three  cases,  the  entrance  foil  scattering  angle  was  A9  =  0.031 , 0.045  and  0.063,  respectively, 
corresponding  to  carbon  foil  thicknesses  of  5, 10  and  20  mils  for  the  10  MeV  MPPE  beam. 
As  expected,  the  thicker  foil  produces  a  significantly  larger  beam  radius,  but  the  tailoring  is 
good  in  all  three  cases. 
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in.  MULTI-STAGE  CONDITIONING  CELLS 


As  we  have  seen,  passive  IFR  cells  can  be  quite  effective  in  tailoring  the  beam 
emittance  so  as  to  reduce  the  growth  rate  of  the  hose  instability  upon  subsequent  injection 
into  the  atmosphere.  However,  the  possibility  of  instability  growth  remains,  albeit  at  a  lesser 
rate.  Since  the  beam  produced  in  an  accelerator  typically  has  significant  transverse 
fluctuations  which  can  provide  the  seed  for  hose  growth,  it  is  also  necessary  to  center  the 
beam  to  danq>  this  type  of  noise.  In  order  to  do  dris,  it  is  necessary  to  pass  die  beam  through 
a  centering  cell  either  before  or  after  the  IFR  tailoring  cell.  In  this  section,  we  shall  discuss 
various  options  for  centering  the  beam,  and  the  ways  in  which  their  compatibility  with 
tailoring  cells  can  be  optimized.  We  shall  not  study  the  centering  process  itself  (all  of  the 
simulation  studies  in  this  pqier  are  axisymmetric),  but  we  shall  infer  the  requirements  for  and 
consequences  of  the  centering  process.  We  shaU  consider  three  types  of  centering  cells:  (A) 
laser-ionized  IFR  charmels,  (B)  passive  IFR  cells  with  narrow  walls,  Hlled  with  gas  at 
relatively  high  density,  and  (C)  wire  cells. 

A.  Preionized  Channel  Centering  Cell  with  Passive  IFR  Tailoring  Cell 

In  some  of  the  ATA  experiments,  a  KrF  laser  was  used  to  ionize  benzene,  in  order  to 
create  an  IFR  channel  that  was  used  to  center  a  beam.^’*’^®'^*’^*  In  this  scheme,  the  gas 
density  is  chosen  low  enough  (<  1  mtorr)  so  that  further  ionization  of  the  gas  by  the  electron 
bean .  can  be  neglected.  This  beam  centering  technique  was  at  times  employed  inside  the 
accelerator  as  well  as  in  a  separate  centering  cell.  A  passive-IFR  tailoring  cell,  containing 
gas  at  higher  pressure  was  located  after  the  centering  cell.^’^  The  success  of  this  techmque 
depends  on  keeping  the  beam  radius  narrow  at  the  point  of  injection  to  the  IFR  tailoring  cell, 
and  the  microdivergence  at  this  point  matched  to  the  IFR  pinch  force. 

In  Table  IV,  we  summarize  the  conditions  for  several  FRIEZR  simulations  that 
illustrate  transport  through  the  complete  active  IFR^assive  IFR  conditioning  system.  In  all 
cases,  the  laser-IFR  cell  was  120  cm  long,  with  a  pre-ionized  channel  with  radius 
Sei,  *  0.8  cm.  The  ion  chaimel  density  was  assumed  to  have  an  on-axis  density  n^^  and  radius 
gfc,  and  the  line-integrated  charmel  strength  f,,  s  Nj/N,,  =  (»nioa^h)/(^b</®c).  A  foil  separated 
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the  laser-IFR  cell  finmi  a  passive  IFR  cell  90  cm  in  length,  filled  with  air  at  density  S  mtotr. 
These  conditions  were  chosen  as  representative  of  the  MPPE  experiment  on  AT  A,  and  the 
beam  parameters  were  those  chosen  in  Sec.  n.C.  to  represent  the  ATA  beam,  except  that 
Cq  =  0.36  rad-cm.  (In  reality,  this  configuratimi  was  only  used  on  earlier  ATA  e]q>eriments^ 
which  had  beam  energies  of  40  MeV;  an  active  IFR  cell  was  only  used  on  MPPE  in 
conjunction  with  a  focusing  lens  and  vacuum  transport  region^  of  the  sort  described  in  the 
previous  section.)  The  results  are  summarized  in  Table  IV  and  Figs.  14-16. 

TABLE  IV.  Summary  of  FRIEZR  simulations  using  a  120  cm  long  laser-IFR  centering  cell 
followed  by  a  90  cm  long  passive  IFR  tailoring  cell  with  air  at  5  mtorr.  The  radius  is  an 
average  beam  radius  in  the  middle  of  the  IFR  centering  cell  (z  =  60  cm).  The  beam  head, 
body  and  tail  radii  at  the  exit  of  the  IFR  tailoring  cell  (z  =  210  cm)  are  denoted  by  the 
quantities  and  a^(40)  which  ate  taken  at  r  =  4, 20  and  40  nsec,  reflectively. 


Case 

U 

cm 

av2(4) 

cm 

av2(20) 

cm 

a,/2(40) 

cm 

Comments 

3-A 

0.005 

0.032 

1.20 

3.4 

1.93 

1.48 

f^  too  small 

3-B 

0.05 

0.032 

0.58 

3.1 

2.2 

0.68 

Near  optimum  f^ 

3-C 

0.15 

0.032 

0.42 

4.0 

3.1 

1.08 

f^  too  large 

3-D 

0.10 

0.045 

0.39 

3.5 

2.4 

0.67 

=  0.24  rad-cm 

3-E 

0.05 

0.045 

0.58 

3.7 

2.8 

0.92 

Thicker  foil 

3-F 

0.02 

0.032 

0.54 

2.9 

2.08 

0.68 

Near  optimum  f^ 

3-G 

Var. 

0.032 

0.47 

2.9 

1.85 

0.78 

f„  =  0.05-0.01 

("Laser-timed") 

Figures  14  and  IS  show  at  the  end  of  the  tailoring  cell  (z  =  210  cm),  and  the 

emittance  after  passing  through  an  exit  foil,  for  three  cases.  These  figures  illustrate  die  effect 
of  vaiying  the  ion  channel  strength  f^  in  the  centering  cell.  For  ^  =  0.005  (Case  3-A  of  Table 
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rV),  the  ion  channel  is  so  weak  that  the  beam  expands  to  an  undesirably  large  radius 
ai/2  s  1.1  cm  in  the  centering  cell.  It  is  undeipressured  at  injection  to  the  tailoring  cell.  The 
beam  tail  expands  a  bit  more  in  the  tailoring  cell  (due  to  the  focusing  effect  of  the  entrance 
foil);  the  head,  with  its  low  emittance,  expands  only  to  =  3.4  cm  by  the  end  of  the  cell, 
and  the  result  is  a  poorly  tailored  profile  with  emittance  ratio  less  than  2:1  after  passing 
through  the  exit  foil.  In  Case  3-B,  the  laser  channel  strength  f^  =  O.OS  is  well  matched  to  die 
subsequent  IFR  cell.  Unlike  the  previous  case,  the  beam  enters  the  tailoring  cell  with 
suitably  small  radius  (~0.6  cm).  The  tail  radius  at  exit  from  the  tailoring  cell  is  0.68  cm, 
close  to  the  injection  value,  while  the  head  radius  increases  to  3.1  cm.  Figure  IS  shows  the 
emittance  is  well  tailored  at  exit,  with  a  ratio  of  3:1  from  head  to  tail.  In  Case  3-C,  the  laser 
channel  is  so  strong,  at  4  s  0.15,  that  the  beam  emittance  increases  considerably  in  the  laser- 
IFR  cell.  Upon  injection  to  the  tailoring  cell,  the  beam  is  overpressured  and  expands 
considerably  in  the  tailoring  cell,  again  leading  to  poor  emittance  tailoring. 

Figure  16  shows  aj^(r)  at  three  different  locations,  for  the  successful  Case  3-B.  At 
z  =  60  cm  in  the  laser-IFR  cell,  the  radius  is  relatively  constant  except  for  an  expanded  bl^  at 
r  s  18  nsec.  This  perturbation  is  due  to  the  nearly  coherent  coU^>se  of  the  ion  channel  at  this 
point.  This  type  of  ion  dynamics  must  be  monitored;  it  can  sometimes  (but  not  in  this  case) 
degrade  the  beam  tailoring.  At  z  =  180  cm  (60  cm  into  the  tailoring  cell),  it  is  evident  that 
the  beam  head  has  not  had  time  to  expand  sufficiendy,  but  by  the  end  of  the  tailoring  cell  at 
z  =  210  cm,  the  tailoring  profile  is  suitable. 

There  are  other  parameters  available  to  optimize  the  matching  of  the  laser-IFR 
centering  cell  and  the  passive  IFR  tailoring  cell.  In  some  experiments  at  ATA,^*^  the  strength 
of  the  laser  charmel  was  reduced  at  the  end  of  the  laser-IFR  cell,  thus  causing  the  beam  radius 
to  increase  slightly  prior  to  injection  into  the  conditioning  cell.  (This  is  done  by  varying  the 
density  of  benzene,  the  laser  target  gas,  as  a  function  of  z.)  The  diverging  character  of  the 
beam,  da|/dz  >  0  as  the  beam  strikes  the  entry  foil  to  the  IFR  tailoring  cell,  can  be  adjusted  to 
improve  the  tailoring.  In  the  ATA  experiments,  this  technique  was  used  to  substantially 
reduce  hose  growth  when  the  beam  was  subsequendy  propagated  in  dense  gas.^  In  Case  3-G 
of  Table  IV,  the  ion  charmel  strength  was  ramped  down  linearly  from  O.OS  to  0.01  between 
z  s  60  and  120  cm.  For  these  particular  parameters,  the  resulting  tailoring  profile  is  found  to 
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be  similar  to  that  obtained  1^  having  a  constant  s  0.02.  At  higher  tailoring  cell  pressures, 
ramping  down  the  channd  strengdi  prior  to  injectitm  into  the  IFR  tailoring  cell  frequently  is 
found  to  alter  the  radius  profile  substantially. 

B.  TwO'Stage  Passive  IFR  Cells 

Femsler,  et  al.,^^  have  calculated  the  centering  of  a  beam  by  a  narrow  high>density 
passive  IFR  cell.  A  cell  of  this  type  has  been  deployed  as  the  fost  stage  of  a  two-stage 
passive  IFR  cell  in  experiments  with  the  SupolBEX  beam.^^  The  configuration  is  shown 
schematically  in  Fig.  2b.  In  these  experiments,  the  centering  cell  (first  stage)  was  80  cm 
long,  with  waU  radius  a^  =  2.5  cm,  and  gas  pressure  =  20  mtorr.  The  tailoring  cell 
(second  stage)  was  40  cm  long,  with  5  cm  wall  radius  and  gas  pressure  3-10  mtorr.  In  some 
cases,  the  two  sections  were  separated  by  a  thin  foil;  in  others,  an  aperture  was  used  in 
conjtmction  with  differential  pttnq>ing  to  maintain  the  pressure  difference  between  the  two 
stages.  We  have  performed  a  simulation  study  of  this  two-stage  conditioning  scheme,  using 
these  device  parameters  and  SuperlBEX  beam  parameters  as  specified  in  die  last  column  of 
Table  I.  We  report  here  on  three  simulation  cases;  furdier  details  wiU  be  rqioned  elsewhere. 

In  Case  4-A,  the  pressure  in  the  tailoring  cell  was  Pi  =  10  mtorr,  and  there  was  no  foil 
(and  therefore  no  scattering  or  focusing)  between  the  two  stages.  Case  4-B  was  similar 
except  that  P2  =  5  mtorr.  Case  4-C  was  similar  to  Case  4-A,  except  diat  die  two  passive-IFR 
zones  were  separated  by  a  foil;  foil  focusing  was  modeled  using  the  usual  thin  lens 
approximation.^ 

Figure  17  shows  the  beam  radius  at  z  =  50  cm  (near  the  middle  of  the  centering  cell) 
and  120  cm  (at  the  exit  of  the  tailoring  cell),  for  the  foilless  10  mtorr  Case  4-A.  There  is 
some  tailoring  even  in  the  narrow  centering  cell,  and  the  continued  evolution  to  a  tapered 
profile  in  the  tailoring  cell  is  evident.  The  final  radius  taper  at  z  =  120  cm  goes  from  2.5  cm 
in  the  beam  head  to  0.7  cm  at  r  =  27  nsec.  The  flat  profile  of  e^iiiir)  for  t  <  6  ns  is 
indicative  of  a  problem:  the  head  of  the  beam  scriq)es  the  wall  of  the  narrow  centering  cell, 
current  is  lost,  and  the  half-current  radius  is  restricted  to  aj/2  “  >  1.8  cm  when  the 

beam  fills  die  narrow  centering  pipe. 
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Figure  18  coiiq>ares  at  the  exit  of  the  IFR  tailoring  cell  for  the  three  cases.  Case 
4- A  shows  the  most  satisfactorily  tailored  profile.  In  Case  4-B,  with  tailoring  pressure 
p2  =  5  mtorr,  the  beam  was  mismatched  with  too  large  an  emittance  at  injection  to  the 
tailoring  cell.  As  a  result,  the  tail  expands,  and  the  profile  of  at  exit  is  rather  flat, 
tiq>ering  only  a  factor  of  two  firom  head  to  tail.  The  inclusion  of  foil  focusing  (Case  4>C) 
has  disastrous  consequences.  In  the  low  current  portion  of  the  beam  head  (1-3  kA),  the 
effective  focal  length  for  the  tailoring  cell  entrance  foil  is  10  to  30  cm;  this  limits  subsequent 
beam  expansion  and  is  responsible  for  the  large  dip  in  the  radius  seen  at  r  ■>  8  nsec.  In  the 
high  current  beam  body,  the  foil  focusing  effects  are  much  stronger,  with  the  focusing  angle 
Bf  ~  0.2-0.3.  This  part  of  the  beam  is  overfocused,  and  expands  rapidly  after  reaching  its 
initial  focus,  as  is  evident  in  Fig.  18  in  the  region  10  nsec  <  r  <  20  nsec.  Indeed,  the  latter 
half  of  the  beam  scrapes  the  wall  of  the  tailoring  cell,  resulting  in  the  loss  of  3  to  4  kA  (in 
addition  to  the  loss  of  approximately  2  kA  in  the  centering  cell)  from  the  original  peak 
current  of  IS  kA.  This  is  qtparent  in  the  plots  of  ^(r),  at  exit  from  the  tailoring  cell,  shown 
in  Fig.  19  for  Ceases  4-A  and  4-C. 

To  avoid  beam  degradation  resulting  forni  foil  focusing,  a  "foilless"  two-stage  IFR  cell 
was  designed  and  installed  on  SuperlBEX.  This  configuration  used  an  aperture  combined 
with  differential  punq>ing  to  generate  the  desired  density  difference  between  the  two  stages. 
It  produced  some  promising  results  but  often  exhibited  substantial  current  loss,  and  did  not 
usually  generate  well-tailored  pulses.  Simulations  of  this  configuration  were  sensitive  to  the 
assumed  scale  length  for  tlK  pressure  gradient  between  the  two  stages;  a  long  gradient  scale 
length  often  caused  substantial  current  loss  near  the  end  of  the  centering  cell  in  the 
simulations.  The  experimental  results  have  been  described  by  Peyser,  et  al.^^  and  Murphy,  et 
al.^ 

C.  IFR  and  Wire  Cell  Hybrids 

In  earlier  experiments  such  as  those  described  in  Secs,  m  A  and  B,  the  beam  traversed 
the  centering  cell  first,  followed  by  the  IFR  tailoring  cell,  which  was  terminated  with  a  thick 
scattering  foil  to  freeze  in  the  emittance  tailoring.  Since  the  beam  is  subject  to  little  or  no 
centering  force  (particularly  in  the  beam  head)  in  a  passive  IFR  tailoring  cell,  die  beam  may 
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drift  off-axis  during  the  tailoring  stage,  thus  amplifying  any  initial  transverse  velocity 
perturbations.  These  perturbations  subsequently  couple  very  effectively  to  hose  instability. 
Thus,  it  would  be  clearly  advantageous  to  center  the  all  of  the  beam  after  it  has  been  tailored. 

Active  wire  cells  provide  an  effective  means  for  centering  a  beam.  In  such  devices,  a 
thin  wire  carrying  a  current  -  (0.1-0.5)I|,  provides  a  strongly  anharmonic  pinching 
magnetic  field  which  centers  the  beam  and  damps  transverse  perturbations  by  phase 
mixing.  Hie  cell  may  be  operated  in  vacuum  or  at  gas  densities  as  high  as  1  atm.  In  vacuum, 
the  wire  produces  the  same  sort  of  pinch  force  on  an  ultrarelativistic  beam  as  a  narrow 
preionized  ion  channel.  (The  laser-IFR  centering  cell  on  ATA  was  often  referred  to  as  a 
"wire  zone"  for  this  reason.)  In  practice,  it  is  convenient  to  operate  a  wire  cell  in  full  density 
air,  but  the  situation  is  then  complicated  by  the  charge  and  current  distributions  induced  in 
the  air,  and  resistive  hose  instabilities  may  be  excited  in  the  wire  cell  if  the  wire  current  is  too 
low.  Wire  cells  operating  in  fu’*  density  air  were  first  developed  by  Frost^  and  have  been 
deployed  on  IBEX  and  RADLAC  at  Sandia  National  Laboratories  and  on  SuperlBEX  at  the 
Naval  Research  Laboratory. These  devices  are  often  referred  to  as  "gas  B ^  cells." 

In  a  wire  cell,  a  beam  reaches  an  equilibrium  with  transverse  temperature  characterized 
by  a  thermal  spread  Ad  «  (21^/7!^)*^  where  I^  2  17  kA.  In  general,  this  represents 
substantial  heating  of  the  beam,  and  is  incompatible  with  the  requirement  of  low  beam 
emittance  on  entry  to  a  passive  IFR  tailoring  cell  unless  the  wire  current  is  much  smaller  than 
the  peak  beam  current.  It  is  not  possible  for  gas  B^  cells  to  operate  in  the  small  I^  regime 
due  to  hose  instability  growth  in  the  wire  cell  itself.  It  is  possible  in  principle  to  operate  a 
vacuum  wire  cell  with  very  small  I^,  but  the  centering  cell  length  would  have  to  be 
substantially  increased,  and  the  problem  of  beam  drifts  in  the  IFR  tailoring  cell  would  persist. 
Thus,  wire  cells  are  generally  unsuitable  for  use  prior  to  an  IFR  tailoring  cell.  However,  wire 
cells  function  very  effectively  when  located  after  an  IFR  tailoring  cell.  In  addition  to 
centering  the  beam,  the  transverse  heating  of  the  beam  by  the  wire  convens  the  radius  taper 
to  an  emittance  taper;  it  is  then  not  necessary  to  use  a  thick  scattering  foU  for  this  purpose. 
Simulations  of  wire  cells  in  both  vacuum  and  gas  indicate  that  the  passive-IFR/wire  cell 
configuratitm  is  very  elective  in  tailoring,  centering,  and  dancing  transverse  fluctuations  on 
the  beam.  These  results  will  be  discussed  in  more  detail  in  a  future  publication.  A  two-stage 
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conditioning  cell  consisting  of  a  passive  IFR  ceU  operating  at  10  mtorr  followed  by  a  gas 
cell  with  •>  5-9  kA  has  recently  been  installed  on  SuperlBEX.^^’^  This  configuration, 
which  is  shown  schematically  in  Fig.  2c,  q>pears  to  produce  well-tailored  pulses  and  shows 
significant  promise  in  reducing  hose  instability  growth. 
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IV.  MULTI-PULSE IFR  CHANNEL  PHYSICS 


Some  potential  q>plications  of  intense  relativistic  electron  beams  require  a  series  or 
burst  of  pulses  widi  sqiaration  time  r, «  10'^  sec.  These  iq[>plications  include  high  energy 
electron  beam  welding  in  dense  gases^^  and  long-range  propagatitm  in  the  atmo^ere.^’^ 
Questions  arise  as  to  possible  effects  which  could  make  a  passive  IFR  cell  unsuitable  for 
multi-pulse  applications.  For  example,  does  die  first  pulse  leave  behind  so  much  plasma  duu 
subsequent  pulses  ate  adversely  affected?  Do  die  eariy  pulses  bore  a  hole  in  die  neutral  gas, 
so  that  subsequent  pulses  produce  a  weaker  ion  channel?  This  section  discusses  the  {diysics 
of  a  passive  IFR  cell  in  multi-pulse  operation  and  addresses  the  fundamental  question  of 
whedwr  the  channel  parameters  will  vary  substantially  from  pulse  to  pulse. 

IFR  channel  physics  for  mult^le  beam  pulses  can  be  thought  of  as  proceeding  in  three 
stages.  The  first  stage  covers  the  period  ^  <  fp  (typically  tens  of  nsec)  during  which  the 
electron  beam  is  present.  During  this  stage,  the  electron  beam  creates  a  plasma  channel  at  a 
rate  Vj  «  10®  P  where  the  pressure  P  is  in  mtorr.^*^  Plasma  electrons  are  quickly  accelerated 
to  energies  of  >100  keV  by  the  space  charge  field  of  the  beam  and  are  radially  expelled, 
leaving  behind  an  ion  chaimel  to  pinch  the  beam.  This  field  also  causes  the  ion  chaiuiel  to 
collapse  toward  the  axis  on  a  time  scale  givoi  by  Eq.  (14).  Typically,  <  Tp,  although 
the  collapse  time  for  an  individual  ion  depends  on  the  time  history  of  h,  and  a^,  after  the  ion 
was  created.  Because  of  the  distribution  of  ion  creation  times,  ion  channel  coUqise  heats  the 
ions  to  a  temperature  of  the  order  of  the  electrostatic  potential  well.  These  ions  may  abo 
produce  significant  inq>act  iotuzation.^  At  the  end  of  the  first  pulse,  the  typical  gas  density 
is  several  orders  of  magnitude  larger  than  the  beam  or  ion  chaimel  densities,  so  the  plasma  is 
still  weakly  ionized.  At  these  densities,  the  plasma  electrons  have  a  mean  free  path  of  1- 
10  m,  so  they  contribute  little  secondary  ionization  before  they  strike  the  wall. 

The  second  stage  occurs  after  the  beam  pulse  has  ended.  A  hot,  positive  ion  colunm 
remains  behind.  In  the  absence  of  collisions  and  neutralizing  electrons,  the  ion  column  radius 
will  typically  expand  to  the  wall  radius  within  100  nsec.  This  effect  has  been  seen  in 
FRIEZR  simulations.^  Charge  exchange  appears  to  be  the  dominant  atomic  collision 
process  in  a  multi-pulse  IFR  cell;  in  some  cases,  the  mean  free  path  for  this  process  may  be 
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comparable  with  the  beam  radius.^^  However,  the  fast  neutrals  which  it  produces  tend  to  be 
lost  to  the  walls  without  significantly  heating  the  bulk  of  the  neutral  gas.  Some  plasma 
electrons  may  be  available  to  neutralize  the  expanding  ion  column,  eq)ecially  if  a  reservoir  of 
available  plasma  exists  at  the  IFR  cell  wall.  However,  FRIEZR  simulations^  of  ion  chaimel 
motion  in  the  Recirculating  Linear  Accelerator,^*^^  which  has  beam  and  channel  densities 
comparable  to  those  in  ATA,  indicate  that  the  ion  column  can  be  expected  to  expand  to  a 
radius  comparable  with  the  tube  radius  a^  before  being  neutralized.  Ibis  occurs  because  the 
hot  neutralizing  electrons  which  are  drawn  toward  the  ion  channel  have  a  line  density 
which  cannot  exceed  the  ion  channel  line  density  Nj  and  are  spread  more  or  less  uniformly 
inside  the  cell.  Thus,  their  space  charge  field  does  not  ppreciably  slow  the  expansion  of  the 
ion  column  until  the  ion  chaimel  radius  is  conparable  to  a^. 

The  third  stage  extends  from  the  time  for  the  ions  to  expand  to  the  wall  (~100  nsec) 
until  the  beginning  of  the  next  pulse  (~1  msec).  The  above  discussion  suggests  that  by  the 
time  the  next  pulse  arrives,  most  of  the  ions  will  have  either  been  lost  or  will  be  spread 
throughout  the  IFR  cell.  If  energy  is  transferred  to  neutrals  via  charge  exchange,  those 
particles  will  presumably  be  lost  to  the  wall.^^  Very  little  hole  boring  will  have  taken  place 
because  all  of  the  relevant  mean  free  paths  are  much  larger  than  the  beam  radius.  Thus,  die 
next  pulse  should  enter  a  cell  whose  neutral  density  is  essentially  unchanged  from  that 
present  during  the  first  pulse. 

In  summary,  if  an  IFR  cell  is  used  to  tailor  a  series  of  beam  pulses  separated  by 
~1  msec,  the  first  pulse  will  create  an  IFR  plasma  column  but  will  not  bore  a  density  channel. 
Substantial  plasma  may  be  produced,  especially  by  ejected  plasma  electrons  striking  the 
walls.  However,  this  plasma  should  not  neutralize  the  ion  column  until  its  radius  is  a 
substantial  fraction  of  the  wall  radius.  The  second  pulse  may  encounter  a  difruse  plasma,  but 
this  should  not  affect  the  beam  provided  the  plasma  density  remains  well  below  that  of  the 
channel  formed  by  the  second  pulse.  Thus,  we  expect  subsequent  pulses  to  be  tailored  in  a 
similar  manner  to  the  first  pulse. 

Beam  radius  da*«  taken  from  the  Multi-Pulse  Propagation  Experiment  on  ATA  suj^rt 
this  view.^  A  series  of  five  pulses  with  t,  «  2-5  msec  was  injected  into  a  passive  IFR  cell, 
and  the  beam  radius  profile  for  each  pulse  was  measured  after  it  exited  the  IFR  cell 


using  liglit  emitted  from  a  quartz  foil.  The  profiles  showed  little  or  no  variation  from  pulse  to 
pulse.^ 
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V.  CONCLUSIONS 


Passive  IFR  cells  offer  a  powerful  technique  for  shiq)ing  the  radius  and  emittance 
profile  of  an  electron  beam  prior  to  injection  into  a  dense  gas.  Even  with  the  simplest 
configuration,  these  profiles  can  be  readily  adjusted  by  varying  the  pressure  in  die  IFR  cell, 
the  thicknesses  of  the  IFR  entrance  and  exit  foils  and  the  length  and  radius  of  the  cell. 
Additional  control  may  be  obtained  by  using  a  focusing  lens  ahead  of  the  IFR  cell.  Passive 
IFR  tailoring  cells  can  also  be  used  in  conjunction  widi  other  techniques  designed  to  center 
the  beam  or  damp  initial  fluctuations.  All  of  these  configurations  have  been  simulated  with 
the  FRIEZR  beam  simulation  code  and  have  been  successfully  deployed  experimentally. 

ACKNOWLEDGMENTS 

Many  colleagues  in  the  charged  particle  beam  community  have  contributed  to  this 
work.  We  wish  to  thank  the  ATA  group  at  die  Lawrence  Livermore  National  Laboratory, 
including  Frank  Chambers,  William  Fawley,  Kenneth  Struve,  Ted  Orxechowski,  Simon  Yu, 
William  Sharp  and  John  Weir,  for  numerous  discussions  and  for  sharing  data  prior  to 
publication.  We  also  thank  the  NRL  SuperlBEX  experimental  group,  including  Robert 
Meger,  John  Antoniades,  Thomas  Peyser,  Donald  Murphy,  Robert  Pechacek,  Matthew  Myers 
and  Jonathan  Santos,  with  whom  we  have  maintained  a  close  collaboration.  We  also 
acknowledge  numerous  conversations  with  Charles  Frost,  who  pioneered  many  of  the  beam 
conditioning  techniques  now  used  on  SuperlBEX,  and  with  Dale  Welch,  Khanh  Nguyen, 
Leon  Feinstein,  John  Wagner,  Brendan  Godfrey  and  Steven  Swanekamp  who  have  carried 
out  theoretical  or  simulation  studies  similar  to  our  own.  Finally,  we  acknowledge  the 
contributions  of  Jonathan  Krall  to  the  FRIEZR  simulation  code  and  A.  Wahab  Ali  to  our 
understanding  of  the  atomic  physics  processes  discussed  in  Section  IV. 

This  worir  was  stqrported  by  the  Defense  Advanced  Research  Projects  Agency,  ARPA 
Order  No.  7781,  monitored  by  the  Naval  Surface  Warfare  Center  and  was  also  supported  by 
the  Office  of  Naval  Research. 


28 


APPENDIX:  FRIEZR  SIMULATION  CODE 


A.  General  Description 

The  simulations  reported  in  this  paper  were  performed  using  the  FRIE23(  code/^ 
FRIEZR  is  a  particle  simulation  code  in  which  particles  are  advanced  in  time  t  with  a  fully- 
relativisdc  particle  pusher.  Fields  are  derived  from  solutions  to  the  full  Maxwell’s  equations. 
Ultrareladvistic  and  paraxial  approximations,  which  are  often  employed  in  beam  codes,  are 
not  used  in  FRIEZR.  However,  the  axial  variable  is  r  s  =  t  -  zjc,  rather  than  z,  and  die  axial 
grid  is  divided  into  slices  with  uniform  width  Ar.  Relativistic  beam  particles  have  nearly 
constant  r  and  thus  remain  within  the  same  beam  slice  for  many  time  steps.  The  code  is 
organized  so  that  only  particles  from  one  r-slice  are  in  core  at  any  one  time;  careful 
bookkeeping  and  I/O  management  are  thus  required  to  keep  track  of  particles  slipping  from 
one  slice  to  the  next.  The  field  equations  are  solved  as  functions  of  r,  r  and  t,  but  beam  and 
plasma  particles  are  pushed  in  Cartesian  coordinates.  A  Courant  condition  exists  which 
requires  the  time  step  At  and  axial  grid  step  Ar  to  satisfy  At  ^  Ar/2.  Typical  conditioning 
cell  simulations  enqiloy  100-400  beam  slices  with  500-1000  beam  particles  per  slice.  The 
original  code  is  described  in  more  detail  by  Krall,  et  al.^^  This  ^>pendix  describes  features 
which  have  been  added  specifically  to  treat  IFR  conditioning  cells.  Note  that  Ref.  47  uses 
S  s:  cr  as  the  axial  variable. 

The  dynamics  of  plasma  electrons  and  ions  occur  on  very  disparate  time  scales.  For 
high  current  electron  beams,  plasma  electrons  are  quickly  accelerated  radially  with  velocities 
which  q>proach  c  and  are  thus  expelled  from  the  beam  on  a  time  scale  ~a|/c  <0.1  nsec.  Ions 
collapse  inward  on  a  time  scale  that  is  slower  by  a  factor  (mjAn^)^^.  FRIEZR  provides  for 
three  options  for  treating  plasma  particle  motion:  (l)Use  particle  simulation  to  compute 
electron  and  ion  flow  (often  very  expensive).  (2)Freeze  the  ions  and  simulate  the  electron 
motion  (appropriate  for  short  time  scales).  (3)Assome  plasma  electrons  are  instantaneously 
ejected  and  use  simulation  particles  to  follow  the  ion  motion.  The  latter  approach  has  been 
used  in  the  simulations  reported  in  this  paper.  This  assumption  makes  it  much  more 
economical  to  treat  the  long  beams  which  are  used  in  present  experiments  but  does  not  treat 
inductive  erosi<Mi.^^ 
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Creation  of  plasma  by  beam  impact  ionization  is  an  essential  feature  of  passive  IFR 
conditioning  cells.  In  the  studies  repotted  here,  FRIEZR  assumes  an  ionization  cross  secticm 
0|  =  10'^*  cm^  for  relativistic  electrons  impacting  neutral  gas.  A  fixed  number  N,p  of 
simulation  plasma  electrons  and/or  ions  are  created  in  each  axial  grid  cell  at  each  time  step. 
The  plasma  particles  are  created  at  the  instantaneous  locations  of  simulation  beam 
particles  which  have  been  selected  at  random.  The  charge  of  the  secondary  particles  is 
weighted  relative  to  beam  particles  by 

H  -  o.n  N  ^6t/cN  ,  (A-1) 

8  X  g  sb  sb  sp 

where  5t  is  the  simulation  time  step,  n^  is  the  gas  density,  W,),  =  Ib(T)/Ibo  ^  weight 
associated  with  simulation  beam  particles  in  a  particular  r-slice,  and  N,),  is  the  number  of 
simulation  beam  particles  in  an  axial  grid  cell.  Beam-generated  ionization  causes  the  plasma 
density  to  increase  with  r,  causing  the  beam  head  to  be  less  strongly  pinched  than  the  body. 

The  original  FRIEZR  code^^  treated  propagation  in  a  preionized  plasma  channel  such 
as  that  produced  by  laser  ionization  of  a  low  density  target  gas.  Such  channels  were 
extensively  used  on  FRIE2^  generates  preionized  plasma  channels  by 

bringing  in  new  plasma  particles  into  the  r  ~  0  beam  slice  at  each  time  step.  The  laser- 
produced  channel  density  may  also  be  ramped  up  experimentally  during  the  electron  beam 
pulse  by  turning  the  laser  on  during  the  beam  pulse.  This  "laser  timing"  technique  was 
extensively  used  on  the  AT  A  accelerator  to  fine-tune  the  beam  pulse. Since  laser 
timing  creates  a  chaimel  whose  density  increases  during  the  pulse,  the  resulting  beam  radius 
profile  is  similar  to  that  produced  by  beam  impact  ionization.  FRIEZR  treats  laser  timing  in 
the  same  manner  as  beam-generated  ionization;  however,  electron-ion  pairs  are  created  from 
a  "phantom"  set  of  particle  positions  chosen  to  reproduce  the  radial  profile  of  laser-produced 
ionization. 

Although  z  is  not  an  independent  variable  in  the  FRIEZR  code,  most  time-resolved 
experimental  diagnostics  are  measured  at  a  fixed  z-location.  In  order  to  facilitate 
experimental  comparisons,  FRIEZR  dumps  data  for  each  r-slice  at  specifred  "diagnostics 
stations"  at  fixed  z  locations,  using  the  variable  change  z  =  ct  -  r.  The  primary  quantity  of 
interest  in  beam  conditioning  studies  is  the  half>current  radius  a|^(r),  defined  as  the  radius 
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enclosing  50%  of  the  simulation  particles  ronaining  in  Oie  system  at  a  given  z-location.  The 
ims  radius  agg,  is  less  useful  because  it  may  be  strongly  influenced  by  a  "halo"  of  particles  at 
a  large  radius.  However,  FRIEZR  calculates  the  emittance  using  die  formal  definititm  of 
Lee  and  Cooper^^  which  is  based  on  the  rms  radius.  The  code  also  calculates  a  variety  of 
plasma  and  electromagnetic  fidd  quantities. 

FRIEZR  allows  the  wall  radius,  gas  density,  preionized  chaimel  density  and  axial 
magnetic  field  to  be  varied  in  a  piecewise-linear  fashion  with  z.  This  capability  makes  it  easy 
to  treat  complicated  configurations  including  hybrid  conditioning  cells  such  as  those 
described  in  Section  m.  hi  addition,  FRIEZR  can  include  discrete  beamline  elements  such  as 
foUs  and  thin  magnetic  lenses.  These  may  be  located  at  any  axial  position  z  s  z^.  Magnetic 
lenses  are  treated  by  specifying  a  focal  length  As  it  passes  through  the  lens,  each 
simulation  particle  is  given  an  impulse  Spg/Po  ~  ^Py/Po  ~  ~y/^L  ^here  x  and  y  are 

the  particle  coordinates  at  z  =  Z|^  and  p^  =  is  the  nominal  beam  momentum  at  injection. 

Magnetic  lenses  were  extensively  used  to  fine-'tune  the  beam  uqier  on  ATA.^ 

Foils  are  widely  used  to  separate  different  beamline  sections.  FRIEZR  treats  scattering 
by  a  foil  by  specifying  a  characteristic  scattering  angle  based  on  the  foil  composition, 
thickness  and  nominal  beam  energy.  For  example,  a  titanium  foil  of  thickness  Tf  (in  mils) 
produces  a  scattering  angle  The  impulse  given  individual  simulation  particles 

is  chosen  at  random  from  a  gaussian  distribution  with  width 

In  addition,  a  conducting  foil  produces  a  focusing  impulse  similar  to  diat  produced  by  a 
thin  lens;  this  is  because  the  radial  electric  field  E,  produced  by  the  beam  is  shorted  out  by 
the  foil.  This  process  has  been  modeled  by  Adler,^^  Humphries  and  Ekdahl^^  and  Femsler,  et 
al.^^  FRIEZR  employs  a  variation  of  Adler’s  thin  lens  approximation  to  treat  foil  focusing. 
The  impulse  has  the  form  given  by^ 

&  -  -6p  /ymc  *  F  (r/a  )  (r/a  )  I  /7I  ,  (A-2) 

f  r  ZD  D  D  o 

The  function  Ff(r)  is  typically  1-3  and  varies  weakly  with  r  and  a^/a^  and  is  slightly  different 
for  different  assumed  current  density  profiles.  The  effective  focal  length  of  this  thin  "lens"  is 
approximately  For  high  current,  moderate  energy  beams  such  as  those  produced  by 
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SuperlBEX,  the  focusing  angle  Bf  can  be  quite  steep,  and  the  focal  length  may  be  only  a  few 
centimeters.  This  "overfocusing"  effect  may  degrade  the  beam  significantly.^** 

A  three-dimensional  version  of  FRIEZR,  called  ELB  A,**^2*“  is  capable  of  treating  a 
variety  of  effects  such  as  ion  hose^^’^*^^  and  electron  hose*  instabilities  which  are  neglected 
by  FRIEZR.  Preliminary  results  suggest  that  these  effects  are  not  important  in  passive  IFR 
cells. 
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HGURE  CAPTIONS 


FIG.  1.  Schematic  diagrams  of  IFR  tailoring  cell  configurations.  Sinq>le  passive  IFR  cell  (a) 
uses  a  thin  entrance  foil  (FI)  and  a  thicker  exit  foil  (F2).  the  region  between  the  two  foils 
contains  a  low  density  gas.  In  configuration  (b),  the  IFR  cell  is  preceded  by  a  vacuum 
transport  region  with  a  focusing  lens  to  control  the  rate  of  beam  expansion  at  the  entrance  foil 
(FI). 

FIG.  2.  Schematic  diagrams  of  IFR  tailoring  cells  combined  with  various  centering  cells. 
The  active  IFR  /  passive  IFR  configuration  (a)  employs  a  centering  cell  with  a  plasma 
channel  (CH)  created  by  a  laser  or  low  energy  electron  beam.  The  two-stage  passive  IFR  cell 
(b)  uses  a  passive  IFR  cell  with  a  higher  gas  density  and  smaller  wall  radius  to  center  the 
beam  prior  to  injection  into  the  tailoring  cell.  The  IFR  /  wire  cell  hybrid  configuration  (c) 
uses  a  current-canying  wire  (W)  to  center  the  beam  after  it  has  been  tailored.  The  foils  (Fl- 
F3)  may  be  relatively  thin.  The  wire  cell  may  be  operated  in  vacuum  or  in  full  density  air. 

FIG.  3.  Normalized  emittance  e„(r)  at  z  =  90  cm  for  Cases  1-A,  1-B,  1-C,  1-D  and  1-E. 
Emittance  for  Cases  1-A,  1-B  and  1-D  remain  near  their  injected  values  except  for  cooling 
due  to  particle  loss  in  the  beam  head  seen  in  Case  1-D.  Cases  1-C  and  1-E  both  show 
emittance  growth  in  beam  body  due  to  oveipinching  by  the  strong  ion  channel. 

FIG.  4.  Comparison  of  beam  radius  from  fixed  ion  FRIEZR  simulation  Case  1-A  with  the 
corresponding  analytical  solution  from  Eqs.  (5)  and  (6). 

FIG.  5.  Beam  half-current  radius  aj^(T)  at  z  =  30, 90  and  150  cm  for  Case  1-B.  The  beam 
head  expands  freely  while  the  body  is  nearly  in  equilibrium. 

FIG.  6.  Beam  expansion  rate  da/dz  versus  r  at  z  =  30, 90  and  150  cm  for  the  case  shown  in 
Fig.  5. 
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FIO.  7.  Plot  of  beam  half-current  and  rms  radii  and  a^,  versus  r  at  z  s  90  cm  from 
FRIEZR  simulation  Case  1-B.  IFR  cell  pressure  is  S  mtorr. 

FIG.  8.  Plot  of  ai/2(r)  and  a^g(r)  at  z  s  90  cm  for  the  20  mtoir  IFR  cell  Case  1-C.  Beam 
radius  ccmtracts  more  quickly  in  die  beam  head  and  is  essentially  cmstam  for  r  >  20  nsec. 

FIG.  9.  Evolutitm  of  beam  radius  aj^(z)  at  r  =  4, 20  and  40  nsec  from  simulatitm  Case  2-D. 
Nominal  focal  length  for  the  upstream  lens  is  1 14  cm.  The  slice  in  the  beam  head 
(r  s  4  nsec)  reaches  its  minimum  radius  near  the  geometric  focal  point,  while  later  slices 
focus  further  downstream  due  to  space  charge  effects.  This  results  in  a  relatively  large 
minimum  beam  radius  at  the  IFR  exit  foil.  Note  that  this  figure  was  reconstructed  from 
ai/2(r)  data  dumped  at  a  small  number  of  z-locations,  so  the  gridding  in  the  code 
(cAt  s  1.5  cm)  is  much  smaller  than  that  indicated  in  the  figure. 

FIG.  10.  Beam  radius  for  the  same  three  slices  for  a  simulation  with  f|_^  reduced  to 
90  cm  (Case  2-H).  This  causes  the  beam  body  to  focus  near  the  IFR  entrance  foil,  resulting 
in  a  smaller  minimum  radius  at  die  IFR  exit  foil. 

FIG.  11.  Beam  radius  profile  tiipir)  at  the  IFR  exit  foil  for  upstream  focal  length  f|^  s  70, 90 
and  1 14  cm  (Cases  2-L,  2-H  and  2-D).  Optimal  tailoring  is  exhibited  by  the  f^^  =  90  cm  case. 

FIG.  12.  Emittance  profile  eg(r)  at  the  IFR  exit  foil  for  the  three  cases  shown  in  Fig.  11. 

FIG.  13.  Effect  of  the  IFR  entrance  foil  scattering  angle  on  the  beam  radius  profile  at 
the  IFR  exit  foil.  Scattering  angle  =  0.031, 0.045  and  0.063  for  the  three  simulations 
(Cases  2-G,  2-H  and  2-J). 

FIG.  14.  Beam  radius  9iip{r)  at  the  end  of  the  tailoring  cell  for  three  different  values  of  die 
laser-IFR  cell  channel  strength  f^.  Simulatimi  cases  are  3-A,  3-B  and  3-C.  Best  results  are 
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obtained  with  the  inteimediate  value  s  0.05  (Case  3-B).  Both  strong  and  weak  channel 
cases  result  in  a  laiger  final  radius  and  a  less  desirable  tailoring  profile. 

FIG.  IS.  Emittance  profile  Cn(r)  at  the  end  of  the  tailoring  cell  for  the  three  simulations 
shown  in  Fig.  14. 

FIG.  16.  Plot  of  beam  radius  a|^(r)  at  three  different  locations  from  a  simulation  of  a  laser- 
IFR  centering  cell  followed  by  a  passive  IFR  tailoring  cell  (Case  3-B).  The  z  s  60  cm  curve 
is  in  the  middle  of  the  centering  cell,  while  the  other  two  locations  are  at  possible  locations 
for  the  IFR  exit  foil. 

FIG.  17.  Simulation  of  a  SuperlBEX  "foilless”  two-stage  IFR  cell  with  =  20  mtorr  and 
P2  =  10  mtorr  (Case  4-A).  The  z  =  50  cm  location  is  near  the  middle  of  the  centering  cell, 
while  z  s  120  cm  is  at  the  end  of  the  tailoring  cell.  The  final  radius  profile  a2y2(T)  exhibits 
fairly  good  tailoring,  but  has  a  substantial  region  in  the  beam  head  where  it  is  relatively  flat. 

FIG.  18.  Coii9>arison  of  the  beam  radius  profile  at  the  tailoring  cell  exit  for  two-stage  IFR 
simulations  4-A,  4-B  and  4-C.  Case  4-B  has  P2  reduced  to  5  mtorr,  resulting  in  a  relatively 
flat  tailoring  profile  and  a  larger  minimum  radius.  Case  4-C  is  a  repeat  of  the  P2  =  10  mtorr 
Case  (4-A)  but  includes  die  effects  of  foil  focusing.  The  inclusion  of  foil  focusing  degrades 
the  beam  radius  profile  significandy. 

FIG.  19.  Comparison  of  the  transported  current  ^(r)  at  the  tailoring  cell  exit  for  simulations 
with  and  without  foil  focusing  (Cases  4-C  and  4-A).  The  current  loss  seen  in  presence  of  die 
foil  is  due  to  an  overfocusing  effect  which  causes  the  beam  to  expand  riqiidly  in  the  taOoring 
cell  following  its  initial  colli^ise  after  passing  through  the  foil. 
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(c)  IFR/Wire  Cell  Hybrid 
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